Atom-atom correlations in time-of-flight imaging of ultra-cold bosons in optical 

lattices 
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We study the spatial correlations of strongly interacting bosons in a ground state, confined in 
two-dimensional square and three-dimensional cubic lattice. Using combined Bogoliubov method 
and the quantum rotor approach, we map the Hamiltonian of strongly interacting bosons onto 
U(l) phase action in order to calculate the atom-atom correlations decay along the principal axis 
and a diagonal of the lattice plane direction as a function of distance. Lower tunneling rates lead 
to quicker decays of the correlations, which character becomes exponential. Finally, correlation 
functions allow us to calculate quantities that are directly bound to experimental outcomes, namely 
time-of-flight absorption images and resulting visibility. Our results contain all the characteristic 
features present in experimental data (transition from Mott insulating blob to superfluid peaks, 
etc.), which emphasizes the usability of the proposed approach. 
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I. INTRODUCTION 

During the last years enormous progress was made in 
the experimental study of cold atoms in optical lattices 
[![. The great advantage of optical lattices as analog 
simulators of strongly correlated Hamiltonians lies in the 
ability of optical lattices to accurately implement lat- 
tice models without impurities or defects. Furthermore, 
ultra-cold atoms confined in optical lattice structure pro- 
vide a very clean experimental realization of a strongly 
correlated many-body problem Strong correlation 
effects, which imply enhanced quantum fluctuations, are 
playing an increasingly important role in recent exper- 
iments on dilute quantum gases [3]. To underline the 
importance of correlations and fluctuations involved the 
studies of correlation functions are called for. One strat- 
egy for boosting the importance of correlations and fluc- 
tuations involves the study of the atomic correlators as 
a function of control of coupling parameters. The atom- 
atom correlation function is defined by: 



C(r,r') = (a(r)a+ (r')> , 



(1) 



where a (r) is a bosonic operator and (...) is statistical 
averaging [4]. This quantity is required for computing 
time-of-flight absorption images, as obtained with ultra- 
cold atoms released from optical lattices Q . We will 
calculate this quantity using the quantum rotor approach 
developed in our previous works 0, § , however now sub- 
stantially supplemented by the implementation of the Bo- 
goliubov method [9|. 

In experiment, one can envisage to measure the mo- 
mentum distribution of the out-going atoms by taking 
a selective time of flight image (TOF). For homoge- 
neous case, time of flight measurements probe the single- 
particle Green's function at equal times i.e. the one-body 
density matrix. In an experiment of time-of-flight imag- 
ing, the cloud of ultra-cold atoms is first suddenly re- 
leased from the harmonic trap. After a time of flight t, 
the position of the atoms is proportional to the momen- 
tum of the atoms in the initial cloud. Finally, an absorp- 



tion image of the expanding cloud of atoms is taken by a 
probing laser. The resulting image provides directly the 
distribution of the momentum space n(k). It is our goal 
of the presented paper to calculate the time-of-flight pat- 
terns using the combined Bogoliubov method and the 
quantum rotor approach and show that our approach 
recreates all the characteristic features that are observed 
in experimental settings. The plan of the paper is as fol- 
lows: in Section II, we introduce the microscopic Bosc- 
Hubbard model relevant for the description of strongly 
interacting bosons in an optical lattice. In the next Sec- 
tion, we perform the evaluation of the atom-atom cor- 
relation function by splitting the bosonic field into its 
amplitude and phase. Furthermore, in Section III, by 
implementing the Bogoliubov method to the amplitude 
part of the bosonic field and combining its outcome with 
the quantum rotor approach, we derive analytically the 
explicit expressions for the atom-atom correlation func- 
tion, time-of-flight absorption images and visibility. In 
Section IV, results of our calculations are plotted: the 
spatial dependence of the correlation function as a func- 
tion of lattice distance R, time-of-flight images and visi- 
bility for various model parameters. Finally, we conclude 
in the Section V. 



II. MODEL HAMILTONIAN 

From a theoretical point of view, description of the 
bosons in optical lattice can be achieved through the def- 
inition of a microscopic Hamiltonian, which can capture 
the main physics of these systems: the Bose Hubbard 
Hamiltonian. Within this model, the bosons move on 
a lattice within a tight-binding scheme and correlation 
is introduced through an on-site repulsive term, since in 
real Bose gases the interaction between atoms cannot be 
neglected in the physical description of the gas. We con- 
sider a second quan tized, bosonic Hubbard Hamiltonian 
in the form (H[IO|: 
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H = -tJ2 [a f (r) a (r') + a+ (r') a (r)] 

(r,r') 



(2) 



The constant t represents nearest neighbors tunneling 
matrix element and is responsible for the dynamical hop- 
ping of bosons from one optical lattice site to another. 
During a jump between two neighboring sites a boson 
gains energy t. The constant U is the strength of the 
on-site repulsive interaction of bosons. Adding a boson 
to already occupied site costs energy U. Furthermore, 
Jt = fj, + yi where /i is a chemical potential controlling 
the average number of bosons. The operators (r) and 
a (r') create and annihilate bosons on sites r and r' of 
a regular two-dimensional (2D) lattice with the nearest 
neighbors hopping denoted by summation over (r, r') . A 
total number of sites is equal to N and the boson number 
operator n (r) = a' (r) a (r) . We use the Hamiltonian in 
Eq. ([2]) to describe a homogeneous (translationally in- 
variant) system, omitting the effect of the external mag- 
netic potential that is usually superimposed on top of 
the optical lattice potential in order to additionally trap 
the atoms. The external potential can be included in the 
Hamiltonian as ^ r e (r) n (r) and would couple to the 
chemical potential term. We discuss such a scenario and 
its consequences in the Section IV. The realization of the 
Bose-Hubbard Hamiltonian using optical lattices has the 
advantage that the interaction matrix element U and the 
tunneling matrix element t can be controlled by adjust- 
ing the intensity of the laser beams. The Hamiltonian 
and its descendants have been widely studied within the 
last years. The phase diagram and ground-state proper- 
ties include the mean- field ansatz [T^, strong coupling 
expansions p^ - [l5j . the quantum rotor approach [lj|, 
methods using the density matrix renormalization group 
DMRG f^Hga. and quantum Monte Carlo QMC simu- 
lations [21R24I. 



III. ATOMIC CORRELATION FUNCTIONS 

A. Transformation to the amplitude-phase 
variables 

The statistical sum of the system defined by Eq. @ 
can be written in a path integral form with use of complex 
fields, a (rr) depending on the "imaginary time" < r < 
ft = 1/ksT. (with T being the temperature) that satisfy 
the periodic condition a (rr) = a(rr + j3): 



Z = / [VaDa]e 



-S\a,a] 



where the action S is equal to 
S [a, a] — I dr 



H(r)+^a(rr) — a (rr 



(3) 



(4) 



Now, we are briefly introducing the quantum rotor ap- 
proach, which has already been employed for the calcula- 
tion of the phase diagram of the cold bosons in optical lat- 
tice [7fl ■ The fourth-order term in the Hamiltonian in Eq. 
@ can be decoupled using the Hubbard-Stratonovich 
transformation with an auxiliary field V (rr) : 



e 2 i-rJ0 

W T, r J ( fdT\-^p±+iV(TT)n(rT)\ 



2tt 



(5) 



The fluctuating "imaginary chemical potential" iV (rr) 
can be written as a sum of static V° (r) and periodic 
function: 



V{tt) = V° (r) + V (rr) , 
where, using Fourier series: 



(6) 



OO 



e=i 



with the Bose-Matsubara frequencies are — 2ir£/f3 
and £ = 0, ±1, ±2, .... Introducing the U(l) phase field 
(j>(rr) via the Josephson-type relation (25j: 



(j) (rr) = V (rr) 



(8) 



with <fi (rr) = d<f> (rr) /dr we can now perform a local 
gauge transformation to new bosonic variables: 



a (rr) = b (rr) e 



;0(rr ) 



(9) 



As a result, the strongly correlated bosonic system is 
transformed into a weakly interacting bosons, submerged 
into the bath of strongly fluctuating gauge potentials 
(which interactions are governed by the high energy scale 
of U) . The order parameter is defined by: 



(a (rr))=<6(rr))Vi 



(10) 



However, one should note that a nonzero value of the 
amplitude (b (rr)) is not sufficient for superfluidity. To 
achieve this, also the phase (rotor) variables, must be- 
come stiff and coherent. The phase order parameter is 
defined by: 



(11) 



where (. . . ) ^ is averaging over phase action to be calcu- 
lated in the next subsection. This reflects the fact that 
all the atoms in the condensate have the same phase and 
form a coherent matter wave. Thus the condensate pos- 
sess a well defined phase associated with the concept of 
so-called spontaneously broken U(l) gauge symmetry. 

After the variable transformations the statistical sum 
becomes: 



Z = [VbVb] [V4>] e 



-S\b,b,d 



(12) 
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with the action: 



S [6, M] = £ J* drb (rr) JU (rr) 

-* E / dT {e~^ (rr) ~ ( r '' r )]&(rr) 6 (r'r) +/i.c.} 

<r,r'> ° 

+ E / dr < s ( rr ) & ( rr )> 6 - p] * ( rr ) b ( rr ) 

. JO 



—i 2 (rr) + i— 0(rr) 
2LT V ' U Vy ' 



(13) 



which will be used as a departure point for obtaining 
bosonic and phase-only actions. 



B. Transformation to the rotor representation for 
phase variables 

The statistical sum can be integrated over the phase 
or bosonic variables with the phase or bosonic action: 

S^[(j>] = -In J [DbVb]e- s i l M, 
S b [6,6] = -In f [T>4>] e -5 [^M 



(14) 



to obtain: 



Z = J [V4>] e- s <>W = J [VbVb] e- s "M 



(15) 



In order to calculate the phase-only action, we use the 
following approximation: 



a (it) =6 (it) e lMrr) « 6 e^ (rr) , 



(16) 



where 6o is static bosonic amplitude, which will be cal- 
culated in the next Section. The phase-only action from 
Eq. (fT4")l can be written explicitly: 







Sm = / dr\J2 



</> 2 (rr) . a ■ . . 
— - — - +i—6{yt) 
2U U Vy ' 



J cos [<f> (rr) — <j) (r'r) 



(17) 



<r,r') 



where J = t\bo\ represents the stiffness for the phase 
field. It is clear that the phase action is non-linear in the 
phase variables <j> (rr) and the statistical sum in Eq. ([15)) 
cannot be calculated exactly. Because of the trigonomet- 
ric nature of the phase variables it is useful to introduce 
a new uni-modular collective field z (rr) = e l ^ rT \ which 
will be treated within the quantum rotor approach by 
making use of the following resolution of unity: 



1=1 dzdzS 



z(rr) 



z (rr) 



-i(p(rr) 



(18) 



The unit length constraint for z (rr) variables is imple- 
mented on average by the formula (see, Ref. |26| ): 



N — z (rr) z (rr) 



= / d\e NX ~ x ^* z(rT)z{rT \ 



(19) 

which introduces a Lagrange multiplier A. Substituting 
Eq. (fT8| into Eq. (fT5|) and integrating by the cumulant 
expansion over the phase variables, the partition function 
reads: 



Z = J [VzVz}d\e NX ~ s * [z ' z] , 



(20) 



In the thermodynamic limit (N —> oo), the integral (|2"0)) 
can be performed exactly by the saddle-point method. 
The quantum rotor action: 

rP 

S z [z, z] = V / drdr' [(\8 rl , - tb 2 ) 5 (r - r') 

+5 rr ,K- 1 (r - r')] z (rr) z (r'r') , (21) 

where K^ 1 (r — r') is a phase correlator, which includes 
dynamic effect of the U(l) phase held: 



Furthermore, the average 



r)-i0(rr') 



(•••>o = 



is taken only over non-interacting quantum rotors: 



dr 



(rr) 



2U 



(22) 



(23) 



(24) 



The Fourier transform of the correlator in Eq. 
zero temperature limit reads: 



K- 1 (a*) = ^-U 



IT 



(25) 



where v (x) = x — [x] — h, and [x] is the floor function, 
which gives the greatest integer less then or equal to 
x. Non-zero value of the order parameter VP b signals a 
bosonic condensation that we identify as superfluid state. 
The phase order parameter in the quantum rotor model 
can be written as: 



i 

1 -^b=nJ2J dr(z(rr)z(r T )) 



(26) 



which fixes value of the phase order parameter ip b ■ The 
Lagrange multiplier saddle-point value "sticks" at criti- 
cality to the value Aq given by 



Ac - J(k. = 0) + K- 1 (cj e=0 ) = 0, 



(27) 



and obeys the Eq. (|2"7| in the whole low temperature 
ordered phase. 
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C. Bogoliubov transformation of the amplitude 
bosonic variables 



In the Bogoliubov approximation the complicated 
many-body quartic Hamiltonian is reduced to a quadratic 
one, which can be diagonalized exactly. In order to cal- 
culate the bosonic action in Eq. (fT4"|) . we follow the Bo- 
goliubov approach by splitting the bosonic operator into 
a Bose condensate macroscopic occupation No = b 2 and 
non-condensed fluctuation part bd (r) @ ■ As a conse- 
quence, the original operator splits into a sum and ac- 
cording to Eq. ([H]) at an alternative representation: 

b (rr) =b + b d (rr) . (28) 

Substituting the Eq. ([25} into the action in Eq. (fI3)). 
and neglecting the terms of the order higher than two in 
bd (r) operators, the action reads: 



iS — So + iSi + S 2 , 



(29) 



where the zero, first and the second order terms contain- 
ing interactions within and between both sub-systems 
are: 



S 



Nf3 



U „ ,2 

-tz - [i+ — \b \ 



b b . 



(30) 



The linear part in the fluctuation operators reads: 



Si=J2j Q dT^-zt-Ti + U\b \ 2 ] fcoMrr)} 

+ Y, dr{\-zt-Tl+U\bo\ 2 ] bold (rr)}, (31) 

and finally, the quadratic part: 

S 2 = \Y. [ dTdr ' [bd(rr), b d (rr')] 
1 „ Jo 



2U\b Q \-A 
U\b \ 2 



U\bo\ 



2U\b Q \ 



bd (rr) 
b d (rr') 



tj^fdr [b d (rr) b d (r'r) + h.c] 



where: 



.4 



M =F 5 (t - r') 



ch 



(32) 



(33) 



Furthermore, the value of bo amplitude can be calculated 
in the saddle point: 



dS N + -, TT\h I 2 n 
= — tz — fi + U \bo\ = 0, 



(34) 
(35) 



dbo 

which results in: 

This implies that the linear term in Eq. (|2H|) is S\ = 
and the bosonic action is given by: 

S\b,b] =S +S 2 . (36) 



D. Correlation function 

Summarizing up the results of the preceding sections, 
the correlation function in Eq. ^ can be written as a 
product of two correlation functions of amplitude and 
rotor fields: 



C(R) = C 2 (R)C fc (R) 



(37) 



where 



C z (R) = C z (rr; r'r) = (z (rr) z (r'r)), 

C b (R) = C b (rr; r'r) = (b (r) b (r')) fe . (38) 



The averagings appearing in Eq. (|38[) are defined by: 
f [VzDz] . . . e~ s ^ 

J [DbVb]...e- Sb ^ 



J [VbVb] e 



-S b \b,b\ 



(39) 



where S z \z, z] and Sb \b, b\ are given in Eqs. (|2 1 1) and 
([56]) . respectively. As a result, the correlation function 
in the real space splits into a product of averages from 
bosonic and phase sectors. The Green's function in the 
bosonic sector reads: 



G b Qsu m ) = b 2 + Gf d (kw m ) 



(40) 



where: 



Gbd (kw n 



2t(£ -e k ) + [/|bo| 2 -iut 
U\bo\ 2 



U\bg\ 2 



El-(iuje) 2 



2t(eo-s k ) + U\bo\-+iu e 



(41) 



and 



£ k = y2i(|-e k ) [2t(|-ek) +2U\b \ 2 ] (42) 

with a dispersion for a simple cubic lattice: 

£k = cos (ak x ) + cos (ak y ) + cos (ak z ) . (43) 

On the other hand, the phase sector leads to the phase 
Green's function: 



G<f,(koj m ) = (z (kw m ) z (kw m )) 

1 



(44) 



Xo-2tb 2 e k + K- 1 (u rn )' 
As a result, the correlation function finally reads: 

C(R) =G (R) [|6 | 2 + G bd (R)l . (45) 
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E. Spatial dependence of the correlation functions 



The correlation function in Eq. ([T]) can be calculated 
along the lattice axis (x or y), or along the diagonal of 
the square lattice: 



R-parallcl — , 0) , 

Rdiagonal — {tICLx j Titty ) , 



(46) 



where a x , and a y are lattice constants in x, and y di- 
rection, respectively and n is an integer. The results for 
two-dimensional square lattice are presented in Figs. [T]i 
and b for the Mott state within the first and the fifth 
Mott lobe, respectively. The correlation function is plot- 
ted in logarithmic scale. As expected, lower tunneling 
rates lead to quicker decays of the correlations. Also, the 
character of the decay becomes exponential as on-site re- 
pulsion U gets stronger (dependence becomes linear in 
logarithmic scale) . Correlations calculated along the lat- 
tice diagonal are weaker than those parallel to the lattice 
axis (see, Fig . Q}i), which is also in agreement with results 
from Ref. i3). 

From this result, we can deduce that the asymptotic 
spatial behavior of the correlation function C (ft) is given 
by the formula: 




R/a 



Figure 1: (Color online) Correlation function C (R) along axis 
direction for various values of jj for two-dimensional lattice 
inside a) the first and b) the fifth Mott lobe; c) compari- 
son of correlations decay for jj = 0.02 and jj = 0.5 along 
axis (circles) and lattice in-plane diagonal (squares) for two- 
dimensional lattice. 



C(R)=exp(-^), (47) 



where £ stands for the correlation length. On the other 
hand, from Eq. ([4"5j) we see that the critical behavior 
(associated with ordering of the phase variables) is gov- 
erned by G<j, (R) [see, Eq. (|44|) ] - a correlation function 
of the quantum spherical model [28] , for which a univer- 
sal quantum-critical (T = 0) properties imply that the 
correlation length as a function of on-site scattering U, 
close to the critical point, should behave as: 



where U c is the critical value of U. 



IV. TIME-OF-FLIGHT PATTERNS 

The phase coherence properties can be estimated from 
time of flight images. When atoms in the superfluid state 
are released from the trap and optical lattice, phase co- 
herence leads to an interference pattern with interference 
peaks arranged in the order reflecting the optical lattice 
symmetry. However, the Mott insulating phase does not 
have long range phase coherence and does not show the 
interference pattern in the time-of-flight images. To iden- 
tify the ordered state, we concentrate on the momentum 
distribution of particle number n (k) , a quantity of ba- 
sic interest that encodes the strong correlations of the 
system: 

n (k) = C ( R ) e ' kR (49) 

R 

This quantity is measured in time-of-flight experiments, 
which are performed by releasing the atomic cloud from 
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the optical trap potential. Due to the fact that the ex- 
pansion is mainly ballistic, the momentum dependence 
of particle numbers represents the interference pattern, 
which typically has two distinct types of behavior: a 
wide maximum signifying the Mott insulating incoher- 
ent state, and a sharp peaks located in the reciprocal 
lattice vectors, which are considered as a signature of su- 
perfluidity in the system. According to the formulas in 
Eqs. ([371)- (gll), we nave: 



n(k) = (4(r)a k (r)) 



mo 



4 ^2 Gbd ( kw ™) + & o4 ^2 G <t> d ( kcj ™) 

m m 

j^pz ^2 G bd(k'a; m )G0 d (k-k';w m /) 
m Q G bd (k) + boG^d (k) 

+^E GM ( k ') G ^( k ~ k ')' ( 5 °) 



where: 

Gbd (k) = 

G 0d (k) = 



2t{e -e k )+n U BE k 
— 1 H — coth — — — 



coth ^ 




40 k 



(51) 



with: 



where w t = Ht/mwo with wo being the size of the on-site 
Wannier function. Therefore, in order to compare the in- 
terference pattern with experiments, we have to calculate 
n (r). However, the experimental absorption pictures are 
two-dimensional projections of the three-dimensional op- 
tical lattice, i.e. they are column-integrated momentum 
distribution over one axis: 



n,j_ (k x ,ky 



dk z \W(k)\ 2 n(k). 



(55) 



Observing that n (k) = n \G x d (k)], where x=cf>,b (see, 
Eq. 1ST))) and making an expansion: 

G x d (k) = g 0x + g x e (k) + 0[e (k)] , (56) 

where e (k) is given by Eq. (pf5|) , we see that the contri- 
bution: 



dk z \W (k)| 2 cosafc z = ^w Q e^^ 2 « (57) 



is negligible for sufficiently large optical lattice depths, es- 
pecially deep in the Mott regime. Thus, we can consider 
an effective two-dimensional system a good approxima- 
tion of the three-dimensional one. The result is presented 
in Fig. [21 In the superfluid phase, the sharp peaks emerge 
denoting long-range phase coherence. In the Mott phase, 
the momentum distribution becomes a broad, featureless 
maximum. 

To asses the interference pattern, a useful quantity is 
often introduced, which is called visibility and measured 
in many experiments [32 34]: 



(58) 



sx 
u ' 



(52) 



We now turn to the description of the interference pat- 
tern observed after release of the atom cloud from the 
optical lattice and a period of free expansion, where the 
phase coherence of the atoms on the optical lattice can 
be directly probed. The density distribution of the ex- 
panding cloud after time t can be represented as follows 



n(r)= Q)V(k)| 2 n(k) 



(53) 



where m is the atomic mass and t is cloud expansion 
time. The quantity W (k) in the Eq. [S3] is the Fourier 
transform of the Wannier function in the lowest Bloch 
band. Typically, the trapping potential is well approxi- 
mated by a harmonic function so that the envelope has 
the Gaussian form: 



7rV 2 



1 / r 2 
exp j 



(54) 



The minimum and maximum intensities n m i n and n max 
are measured at the same distance from the cloud cen- 
ter, which for two-dimensional square lattice are at k = 
(27r/v / 2a, 2n/y/2a) and k = (2n/a, 0), respectively, at 
which the Wannier envelope in Eq. ([S4")l cancels out. The 
resulting dependence of visibility on interaction ratio is 
presented in Fig. [3J In the superfluid phase, its value 
is almost 1 signaling strong coherence peaks in time-of- 
flight patters. Deeper within Mott lobes, the visibility 
is decreasing, while the time-of-flight pattern becomes a 
circularly symmetric blob. However, in Fig. [5] the phase 
transition occurs abruptly and signature of superfluid- 
ity, namely sharp peaks, disappear fast within the Mott 
phase. This is not usually observed in experiments, where 
the remainders of the peaks are visible even for systems 
deep in the Mott phase. It results from the fact that in 
the experimental setups the atoms are trapped not only 
within the optical potential, but additionally within mag- 
netic parabolic trap that is symmetric around the trap 
center. As a result, the system is no longer homogeneous 
and correlation functions become dependent not only on 
the distance, but also on specific location, which makes 
it impossible to calculate them in analytic way anymore. 
The best we can do is to assume that the homogeneous 
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solution works locally, depending on a "local" chemical 
potential 71 (r) = 7i + e (r) (see, Ref. [35|]). Since, the 
time-of-flight patterns are global absorption images of 
atoms released from the traps, we sum them over each 
lattice site [which has a local value of 71 (r)] multiply- 
ing by a local number of atoms. Because, the value of 
71 (r) is dependent on the distance from the trap center 
(decreasing from its maximal value to zero, on the trap 
boundary), for chosen interaction ratio t/U, parts of the 
system are superfiuid, while the other parts can be in 
the Mott regime. As a results, even if the majority of 
the system is in the Mott state, some remnant signatures 
of superfluidity may still be visible. In order to make 
a reliable comparison of our theoretical prediction with 
experiments, we need a translation of the parameters of 
the Bose-Hubbard model t/U into the quantities deter- 
mining the experimental setup, namely Vq/Er (Vq is the 
potential depth and Er is the recoil energy). The band- 
width parameter t is essentially the gain in kinetic energy 
due to nearest neighbor tunneling. In the limit Vq 3> Er 
it can be obtained from the exact result for the width of 
the lowest band in the ID Mathieu-equation [36]: 



7T \ErJ 



(59) 




The relevant interaction parameter U is thus given by an 
integral over the on-site wave function W (x) via: 



Figure 2: (Color online) Simulation of time-of-flight ab- 
sorption images [from. Eq. (|53[) [ for various interactions 
strength t/U: transition from superfiuid (top-left, t/U = 
0.085) to Mott phase (top-right to the bottom-right t/U = 
0.085, 0.07, 0.055, 0.04, 0.025, 0.01, respectively), for n/U = 
0.5. 



U = 



Ai:h 2 a 



jW(x)| 4 



4V2vr 



a [Vo\ 



3/4 



A \ErJ 



(60) 



where, A and a are laser wave length and scattering 
length, respectively (for 87 Rb, a s — 5.45nm and A = 
850nm). It follows, that: 



t 

U 



V2tt 



^\ -2y/V /E B 



(61) 



After averaging of the time-of-flight images over a range 
of chemical potential /x, we can present our results in 
a form that is directly comparable with experimental re- 
sults in Fig. [4] (see, the figure caption). It is clear that all 
the characteristic features have been reproduced. Also, 
because of the better resolution and lack of noise in the 
theoretical results, hints of superfiuid state that is present 
in narrow spheres around the trap center, are visible even 
for deep lattice Vq/Er = 20. 




Figure 3: (Color online) Visibility calculated from simulated 
time-of-flight images as a function of coupling ratio U/t for 
\x/U = 0.5. 
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Figure 4: (Color online) Simulated time-of-flight patterns 
for comparison with experimental time-of-flight absorption 
images (for values of Vo/Er = 7,10,13,14,17,20, around 
150,000 lattice sites and average number of up to 2.5 atoms 
per lattice site in the center of the trap, chosen as in Ref. [J, 
Fig. 2c-h, ). 



V. CONCLUSIONS 

In summary, we have studied the correlations of cold 
atoms loaded in a two dimensional optical lattice. In 
the superfluid regime the validity of the Bogoliubov ap- 
proximation is restricted to the very weakly interacting 
regime (U/t <C 1) and its breakdown as quantum corre- 
lations become important. In general Bogoliubov quasi- 
particle states correspond to solutions of the approximate 
Hamiltonian with a plane-wave character. As quantum 
correlations become important, the exact eigenstates of 
the interacting system do not necessarily have the sim- 
ple plane wave character, especially in the commensurate 
filling situation where for a critical value of U/t the sys- 
tem exhibits the superfuid-Mott insulator transition. To 
address these issues a new theory beyond the simple Bo- 
goliubov approximation was developed that incorporates 
the phase degrees of freedom via the quantum rotor ap- 
proach to describe regimes beyond the very weakly in- 



teracting one. This scenario provided a picture of quasi- 
particles and energy excitations in the strong interaction 
limit, where the transition between the superfluid and 
the Mott state is be driven by phase fluctuations. Taking 
advantage of the macroscopically populated condensate 
state, we have separated the problem into the amplitude 
of the Bose field and the fluctuating phase that was ab- 
sent in the original Bogoliubov problem. Subsequently, 
the functional formulation this formalism was shown to 
be a powerful tool that incorporates properly the interac- 
tion aspects characteristic of the quantum phase dynam- 
ics. This formalism provides a useful framework, where 
the one particle the correlation functions are treated self- 
consistently and permits us to test and simulate Bose- 
Hubbard Hamiltonian with a whole range of phenomena. 
To explore these a diagnostics tools that include spec- 
troscopy of occupation numbers and correlation measure- 
ment is called for. Especially, the time of flight imaging 
is a very powerful technique to probe the quantum gases 
in optical lattices. From the beginning of the ultra-cold 
atom field, it permitted to observe the Bose Einstein con- 
densation. In this work we have calculated time-of-flight 
patterns and the relevant correlation functions on two 
dimensional optical lattice and compared them with ex- 
perimental data. Our theoretical pictures match the ex- 
perimental ones, where one can observe the formation of 
a central and various neighboring peaks, out of feature- 
less Mott state, which subsequently become sharper in 
the superfluid phase as a result of a phase coherence. It 
would be interesting to map out the nature of the Bose 
transition in optical lattice systems by using further diag- 
nostic tools such as Bragg spectroscopy, which reveals the 
whole momentum structure of the single particle correla- 
tion function. A detailed understanding of this quantity 
is a challenging task in the further experimental studies 
of trapped ultra-cold gases. 
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